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Abstract 



Ph| We present a lattice determination of the leptonic decay constant Fd s of the D s -meson using 
Jjn \ its mass and F K as experimental input. Setting the scale with F K = 160 MeV, our final result 
is Fd s = 252(9) MeV. The error contains all uncertainties apart from the quenched approxi- 
mation. Setting the scale with the nucleon mass instead leads to an decrease of 20(1) MeV 
$_i ■ of the leptonic decay constant. 
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1 Introduction 



The D^-meson consists of a c and a s quark. It is stable in QCD and decays by an emission 
of a W + -boson into a lepton and a neutrino. The amplitude of this decay is characterized 
by the decay constant Fd s which is defined by the QCD matrix element 

(0|A M (0)|D s (p))=ip M F Ds (1.1) 

of the axial current = sj^c. 

Given the CKM-matrix element V^ s , Fn e can be measured experimentally by studying 
the branching ratio BR(D S — > \v). The current status of the experimental determination of 
Fd s has been summarized in PQ. Currently, the most precise data come from the ALEPH 
experiment 0, F Ds = (285 ± 19 ± 40) MeV, and from CLEO 0, F Da = (280 ± 17 ± 25 ± 
34) MeV. In the next couple of years CLEO aims at reducing their error to two percent 
precision 

The status of lattice computations for Fd s |5|6|7j was reviewed in |8|9j with a world 
average of F Ds = 255 ± 30 MeV. With QCD sum rules one currently gets F Ds = 235 ± 

24 Mev mam. 

The goal of this work is a computation on the lattice of the weak decay constant Fn s 
to 3 percent accuracy up to the quenched approximation. All the systematic and statistical 
errors will be analyzed and we will also estimate the size of the quenched scale ambiguity. A 
precise quenched calculation of Fj) s together with precise experimental data supplies us with 
a test of quenched lattice QCD. This is of importance for the determination of the B-meson 
decay constant /b from lattice QCD which is an essential input in the measurement of the 
third generation CKM-matrix elements. 

2 Strategy 

The decay constant Fd s is defined in terms of a QCD matrix element. To evaluate this 
matrix element we have to eliminate the bare parameters of the QCD Lagrangian, that 
is the bare gauge coupling go and the bare masses of the relevant quarks, in favour of 
physical observables. One possible hadronic scheme, which has been summarized in ^2] and 
developed in |13ll4ll2j is to use the decay constant F K = 160 MeV to set the scale 1 and 
eliminate the bare masses of the strange and charm quarks in favour of the masses of the K 
and the D s -meson. Neglecting isospin breaking, the quark mass ratio M s /M = 24.4 ± 1.5 
with M = \{M U + Md) is taken from chiral perturbation theory [To] . 

The results in |13pi4|12] are given in terms of the Sommer scale tq which is derived from 
the force between static colour sources. r /a has been computed as a function of the bare 
coupling g to a high precision for a wide range of cutoff values |16ll7|18j . It is only affected 

M.e. to compute the lattice spacing a in physical units as a function of the bare coupling 
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by lattice artifacts of 0(a 2 ) in the quenched theory [18J. The relation to other hadronic 
scales is known [Hj, r F K = 0.415(9), and r m N « 2.6. This illustrates the inconsistency 
of the quenched approximation since the ratio F K /m N deviates by approximately 10% from 
its experimental value. This is the typical size of the quenched scale ambiguity. Setting 
the scale with Fk = 160 MeV corresponds to setting ro = 0.5 fm while win = 938 MeV 
corresponds to ro = 0.55 fm. Below we will estimate the size of the quenched scale ambiguity 
for the weak decay constant of the D s -meson. 

This work aims at a scaling study of _Fb s on the lattice. Thus we keep the physical 
conditions constant and only vary the lattice spacing. We perform numerical simulations in 
O(a) improved lattice QCD using Schrodinger functional boundary conditions |19|20j on a 
L 3 x T space time cylinder. The practical advantages of this approach are discussed in |21j . 
For unexplained notation we refer to [22j. 

In particular we define the meson sources 

= fl6 E^(y)^(z), C = a 6 £c;.(y) 75 C(z), ( 2 - 2 ) 

with flavour indices i ^ j at the x = and the x = T boundary time slices, respectively. 
The correlation functions 

f A (x )=-±(OM*))> /p(«b) = ~»)), h = -^- 6 (0'0) (2.3) 

are then used to compute the decay constant Fd b at finite lattice spacing. Here A^(x) = 
' l l J i( x )lii.l5' l l J j(. x ) 1S the axial current and P(x) = t/} i (x)j5ipj(x) is the pseudoscalar density. 
The leading O(a) cutoff effects in f A are canceled by using a clover improved action 23 
with coefficient c sw and by improving the axial current, 

A^x) = A^x) + c A ad^P(x). (2.4) 

Both coefficients c sw and c A are known non perturbatively in the quenched approximation 
for bare couplings g% = 6/(3 < 1 [21]. The improved correlation function /j[ is then defined 
in analogy to f A . 

On the lattice the axial current receives a multiplicative renormalization by the finite 
factor Z A . The renormalized axial current is thus given by 

(A R ) M = Z A (l + b A (am qtl + am qj )/2) A 1 ^ + 0(a 2 ), (2.5) 

where the O(a) artifacts that are proportional to the bare subtracted quark masses am q are 
canceled by the term proportional to 6a- 

The correlation functions fj± and /p decay exponentially proportional to e~ x ° muB when 
the excitations due to higher states are small enough. Thus the meson mass m^ s , which has 
the experimental value r mn s = 4.988 | 25 4 , can be obtained from the plateau average of the 
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effective pseudoscalar mass 



am PS (x + a/2) = log J^ X °) (2.6) 

f{{x + a) 

The weak decay constant of the D s -meson at finite lattice spacing can then be obtained 
through 2 

F Ds = -2Z A (l + 6 A (am q , + am qj )/2) A (mDsL 3 ) -i/2 e (, -T/2 )mDs 

v h 

x {1 - rfre-* * - f)V- (T - lo)mG } + 0(a 2 ). (2.7) 

Here the factor (mD s i 3 ) _l/2 takes into account the normalization of one particle states. 
The contribution f 1 cancels out the dependence on the meson sources. Because of the 
exponential decay of the correlation function the product in (J2.7)) is expected to exhibit a 
plateau at intermediate times when the contribution r]^ B e~ x ° A of the first excited state and 
the contribution ^ e -( T - x o) m G f rom the ++ glueball both are small. A plateau average can 
then be performed to increase the signal and is understood in (|2.7)1 . Further explanations 
for equation ()2.7j) and details can be found in [2*T] . 



3 Numerical results 
3.1 Parameters 

Our choice of simulation parameters is shown in tabled The bare couplings and the hopping 



p 


^■meas 


L/a 


L/r 


^critical 




K c 




6.0 


380 


16 


2.98 


0.135196 


0.133929 


0.119053 


4.972(22) 


6.1 


301 


24 


3.79 


0.135496 


0.134439 


0.122490 


4.981(23) 


6.2 


251 


24 


3.26 


0.135795 


0.134832 


0.124637 


5.000(25) 


6.45 


289 


32 


3.06 


0.135701 


0.135124 


0.128131 


5.042(29) 



Table 1: Statistics and parameters for our simulations and demonstration of constant physical 
conditions. 



parameters for the strange quark are the same as used in ^2] , and the hopping parameters for 
the charm quark have been found in that work. The non perturbatively defined improvement 
coefficients c sw and ca are taken from [25- &a has been obtained non perturbatively in [26J. 
Since it is difficult to compute (cmp. |2Zj) and since we have to extrapolate their results 
slightly we have used one-loop perturbation theory j2HI as well. The conversion to physical 

2 In |21| the first minus sign is missing. 
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units is done with the Sommer scale 7*0 . From the fit function of ^1 we infer that our lattice 
size L/a — 16 at (5 — 6.0 corresponds to L/r = 2.98. The results of |Ej show that at this 
volume the finite size effects may be neglected. The lattice spacing a is varied from 0.1 fm 
to 0.05 fm. 

Furthermore we choose T/L = 2 in contrast to T/L = 2.5 used for the computation 
of the charm quark mass in J2| to reduce rounding effects, which would be the dominant 
systematic effects on F Ds otherwise. 

The fourth column of table Q shows that the above choice of simulation parameters 
corresponds to approximately constant physical box sizes. Also the mass of the D s -meson 
takes its correct value up to small fluctuations. We have checked that these lead to negligible 
corrections for the decay constant Fd s . 

3.2 Computation of the D s -meson decay constant 

To compute the meson mass i?id s and the decay constant Fd s we calculate the combina- 
tions ()2.6I2.7|) of correlation functions. For all the parameter choices we find plateaus as 
functions of xq. Our task is to find the plateau region such that the effect of excited states 
is negligible compared to our statistical error. 

When dealing with heavy quark propagators on single precision machines there is always 
the danger of roundoff problems which might be the dominant uncontrolled systematic error. 
Since the computation of the decay constant i<b s involves propagators extending from the 
xq = to the xq = T boundary in the correlation functions f\ the discussion of these issues 
becomes even more important as for the computation of masses. From we estimate that 
for T = 2L instead of T = 2.5L taken there the rounding errors for Fd s should be small 
enough. To show this we computed our observables on O(10) configurations at (3 — 6.0 and 
6.45 in single as well as in double precision. This check reveals that the rounding errors for 
Fd b are of the order of one per mill. 

From fits of Fd s (xq) using an a priori chosen plateau range we extract the contributions 
rf^e~ x °^ and rf A e~^ T ~ XoS>m-G of the first excited state and the ++ glueball excitation, respec- 
tively. The theory in [21] predicts relations between the prefactors and the corresponding 
prefactors for the pseudoscalar masses which are roughly 3 fulfilled. 

All the systematic effects that deteriorate the plateau are added. With an a priori chosen 
upper bound 5 for the maximal relative systematic error in a point Fj) a (xo) we define the time 
interval for the plateau region such that for all times in this interval the added systematic 
contribution is smaller than 5. The systematic error for the plateau averaged observable is 
then even smaller. For the decay constant we took 5 = 0.5%. The plateaus defined in this 
way extend from 4ro to 5ro at the four (3 values considered. For the pseudoscalar mass we 
used 5 = 0.3%. 

3 In some cases the fit function is not determined precisely. Then we choose the fit function such that the 
contribution of the excited states is rather overestimated. Note that we do not attempt to measure A or r/^ s 
but rather want to estimate our systematic effects. 
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3.3 Continuum limit 



We perform the plateau averages of (J2.7|) and obtain the decay constant at finite lattice 
spacings as shown in table 121 These data can be extrapolated to the continuum limit, 



p 


aF^ re 


c A =0 


daF^ re /dc A 


roF Ds 


6.0 
6.1 
6.2 
6.45 


0.1144(25) 
0.0979(17) 
0.0863(20) 
0.0633(14) 


0.2375(58) 
0.1699(33) 
0.1280(32) 
0.0664(17) 


0.540(14) 
0.576(13) 
0.598(16) 
0.614(15) 



Table 2: Simulation results for F Ds . 

leading to our main result. Since we employ O(a) improvement, the natural ansatz for the 
continuum extrapolation is linear in a 2 . This extrapolation is shown in figure ^ in the fit 




Figure 1: Continuum extrapolation of Fd s . 

we neglect the point at (3 = 6.0 which is farthest away from the continuum limit. We quote 

r F Ds = 0.638(24) (3.8) 

as our final result. Using ro = 0.5 fm this corresponds to i*b s = 252(9) MeV. Using 1-loop 
perturbation theory for 6a instead of the results of |2E] leads to 

r F Ds = 0.631(24), (3.9) 
5 



corresponding to Fd s = 249(9) MeV. 

3.4 Quenched scale ambiguity 

To estimate the quenched scale ambiguity we consider r i<b s as a function f(z) of z = r mn s 
around the physical value of z. This is possible since in addition to the hopping parameters 
leading to the physical strange and charm quark masses we have computed our observables 
at five more hopping parameter values around the hopping parameter of the charm quark 
for each bare coupling considered. Here we will only use the two hopping parameters closest 
to the charm value and present further results in a separate publication. 

We expand f(z) = f(z ) + (z — zo)f'(zo) + . . . in a Taylor series around z = 4.988. 
f'(zo) is estimated from a linear fit of roF^ as a function of romD s - A 10% increase of 
r corresponds to z — z = 0.5. Our estimate of the corresponding effect on r F Ds in the 
continuum limit is 

0.5f(z ) = 0.008(3). (3.10) 

Converting back to physical units (now using r = 0.55 fm) we estimate that F Ds decreases 
by 20 MeV corresponding to eight percent. 

4 Conclusions 

Our main result is r -Fb s = 0.638(24). Using r = 0.5 fm this corresponds to Fd s = 
252(9) MeV. This is compatible with all the experimental data currently available and with 
most theoretical estimates as well. The precision of our result however matches the goal of 
precision of future experiments, for example CLEO. 

Under a scale shift of 10%, which is typical for quenched spectrum computations Fd s 
decreases by eight percent. Assuming that all the scale ambiguity comes from neglecting 
dynamical sea quark effects this might also give an idea of the unquenching effect for F^ s 
which however has to be probed in an unquenched simulation. 

Since our complete data set consists also of simulations at larger quark masses we will be 
able to look at the functional dependence of the weak decay constant around the charm quark 
mass. Together with a new computation of the decay constant in the static approximation 
we will then interpolate to the _B-meson mass and compute Fb with a better precision than is 
currently available. The corresponding analysis has still to be done and further simulations 
will be necessary so that these exciting results will be presented in another publication. 

This work is part of the ALPHA collaboration research programme. It was supported 
by the European Community under the grant HPRN-CT-2000-00145 Hadrons/Lattice QCD 
and by the Deutsche Forschungsgemeinschaft in the SFB/TR 09 and the Graduiertenkolleg 
GK271. All the production runs were carried out on machines of the APE1000 series at 
DESY-Zeuthen. The checks of the rounding errors in double precision were done with a C- 
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code based on the MILC collaboration's public lattice gauge theory code, see [29]. This code 
ran on the IBM p690 system of the HLRN |3U] anc ^ on the PC cluster at DESY-Zeuthen. 
We thank the staff at the computer centres for their help, and M. Delia Morte, S. Sint, 
R. Sommer, and U. Wolff for useful discussions. 
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